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^D ' Abstract: Allowing for antineutrinos (vx) and sterile neutrinos {ug) to accompany LMA 

,^ ' oscillations, we derive in a model independent way, upper bounds on their components ""^ 

in the solar flux, using the recent data from SNO and SuperKamiokande. Along with the ^""^ 

^H. general case (LMA + z^^; + z^s) we consider the special cases where only Ux or Us are present. (^ ^ 

r~| , We obtain an upper bound on Px which is independent of the I's component. In the no /- — n 

sterile case we obtain upper and lower bounds on /b, the SSM normalization factor. We 
K^ \ also investigate in the general case the common parameter range for fs and the I'x, Vg 

\-{ ' components and find that the upper bound on Vg is hardly sensitive to the Vx component. 

In the absence of Vx we recover the Ug upper bound existing in the literature. We finally 
present a simple x^ analysis of all four cases considered. 
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1. Introduction 



m 

Despite the recent realization ^ that the solar neutrino deficit, acknowledged over three ^~^ 



survival probability Pee IQ ^^'^ sterile neutrino component bounds |l^], |p^,|l^,|16|,|17|. 
In section 2 we examine the consequences from SNO and SK for the joint possibility of 
active antineutrinos and sterile neutrinos in the solar flux and for the limiting cases of 
each of these components alone. So in this respect our results generalize those of ref. 
p5[| . Regarding Ue, all considerations derived for antineutrinos of the other flavours would 
apply, if not for the recent and very strict upper bound on the former from the KamLAND 
experiment [18|. We find that SK data on the neutrino electron scattering total rate 



leads to the exclusion of all active antineutrinos up to 1.17 a and 0.83(T when combined 
with SNO data from salt phase I and II respectively. Up to 2a the percentage of non 
electron antineutrinos I'x in the active neutrino flux is smaller than 64% and 88% when 
SK is combined with data from SNO I and II respectively. Our results are independent 
of standard solar model normalization. Section 3 deals with a simple x^ analysis and in 
section 4 we expound our main conclusions. 
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J^ 



decades ago [S], results mostly from neutrino oscillations through LMA p], it is by no 
means certain whether oscillations are accompanied by the conversion in the sun of electron 
neutrinos into sterile ones, into antineutrinos of other species, or both. This open question 
has obvious implications in a possible time modulation of the solar neutrino flux, an effect L^-' 

for which evidence was found by the Stanford Group y],[|5[, [^, iQ]. In fact, if electron C ^ 

neutrinos produced in solar fusion reactions interact via a sizeable magnetic moment Q ^~^ 

with a time varying solar magnetic field, the result is the production of a time dependent 
component of active i^^ or u-,- or unobserved sterile neutrinos Q in the neutrino flux from 
the sun, reflecting in some way the time variation of the solar field. 

In this article we perform a model independent analysis of the implications from the v J 

SNO salt phase I and II jl^ and SuperKamiokande (SK) results [jll| on the flux of sterile C ^ 

neutrinos and active antineutrinos which may accompany the LMA effect. Being model /" — ^ 

independent, our analysis will mainly focus on solar neutrino data and its implications 
on bounds of sterile neutrino and antineutrino components. Several model independent 
analyses of solar neutrino data have been performed in the past. These concentrated in the 



2. Model Independent Analysis 

We start with the event rate expressions for the charged current (CC) and neutral current 



(NC) reactions for SNO and neutrino electron scattering (ES) for SK, SNO [19| 



R^^ = fsPee (2.1) 

R^^ = fj^p^^ + /^(l - P^^)[sin^a sin^il) + f^ sir?a cos^iIj] (2.2) 

R^^ = f^Pee + /ij(l - Pee)[r sir? a sir?-i\) + f sir?a cos^^)]. (2.3) 

Of all electron neutrinos that are converted, proportional to 1 — Pee, sir?a denotes the C_J 

fraction that is converted into active ones i^x^i^x, (x ^ e) while ip is the angle describing |~[~] 



the h'xji^x components. The sterile neutrino component is therefore proportional to cos'^a. 
Parameter /b denotes the normalization to the standard solar model ^B neutrino flux |20|. 



l-r) 

Quantities r, r are respectively the ratios of the NC neutrino and antineutrino event rates 

to the NC+CC neutrino event rate and f^ is the ratio of the antineutrino deuteron fission v — J 

to neutrino deuteron fission event rate. Specifically C !) 



/ dE,<P{E,) J dE, f dE'^^fjE',, E,) 
O'NC -^ O-NC+CC 



^.r„ ^ a NC+CC 



J dE^^jE,) J dE, J dE'^^fjE',, E,) Q 



(^NC -^ <^NC+CC 

_ j dE,ct){E,)<jNc{Eu) ,^ „. 

Td = = (2.6) 

0"7VC -^ ctnc 



where / is the energy resolution function for SNO [£^] or SK p^ . Owing to its near energy 
independence in this range, the electron neutrino survival probability Pge is factorized out 
of these integrals as in eqs.(2.1)-(2.3). Energy thresholds considered are E^^^ = 5.5, 5 MeV 
for SNO and SK respectively and the rest of the notations is standard. We obtain 

r = 0.150 (0.151 for SK), f = 0.115 ( 0.116 for SK), fa = 0.954. (2.7) 

where these minor differences are mainly the result of the difference in the threshold ener- 
gies, being largely independent of the resolution functions. The data from SNO (phase I 
and II) and SK are summarized in table I. 



CC NC ES 



SNO I 0.275±[J:[JJ| 0.900 ± 0.081 0.382±[J;[j|^ 

SNO II 0.294±[j:[]20 0M6±lfi 0.368±Hio 

SK 0.406±[j:[|}f 



Table I - Ratios of event rates to standard solar model \2(\] event rates at SNO and SK 
(theoretical error not considered). 
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From eqs.(2.1)-(2.3) one can eliminate the angle a and express angle '0 in terms of the 
experimental event rates and other model independent quantities, 

sm %> = —- — = (2.8) 

7(1 — rd) + r — r 

with 

^^ RNC_RCC- (2.9) 

If not for the large uncertainties that are propagated into equation (8) originated from 
the uncertainties in R^^ , R'-''^' , R^^' , this would unambiguously determine the relative 
proportion of non electron antineutrinos I'x in the active non zvg flux. Hence, as will be seen, 
only upper bounds on the fraction of the Ux flux can be derived at the present stage. To 
this end we evaluate the parameter sin'^ip using eq. (2.8) for all values of R , R , R PJ 

within their allowed la ranges for SNO: these are represented by the light shaded areas l~rj 

in fig.l where the SNO data used are those from salt phase II. Hence for each chosen ^ — n 

value of sin'^ip the allowed values of the three reduced rates lie within each shaded area. 
If the SNO experiment alone is considered, it is seen that all possible values of sin^il) in 
the range < sinPjp < 1 can be obtained. We note that as the i>x component decreases ^""^ 

and eventually vanishes {sin'^'ip -^ 1), the factor multiplying r in eq.(2.3) increases while L nJ 

the one multiplying f approaches zero. Owing to the relative difference between r and (^ ^ 

f [eq.(2.7)], this implies a slight inclination into larger values of R^'^ for a decreasing u^ 
component. The same effect, although much less significant, because of the much smaller 
difference between fd and unity, is also present in the {sin'^'tp, R ) area and does not exist 
inR^'^ [seeeq. (2.1)]. ' ' 

The data from the SK experiment with R^^ restricted to its SK la range (see table C ^ 

I) are also used to evaluate sin^il). In fig.l the dark shaded area which is part of the total (^ ^ 

{sin^ijj, R^'^) one represents the parameter range allowed jointly by SNO II and SK. This 
is enlarged in fig. 2 with a magnified horizontal scale. The result is a lower bound on sm^V 

SNO n sin^^ > 0.5 at 1.49cj(87% CL) (2.10) 

which is in fact an upper bound (0.5) on the fraction of non-electron antineutrinos v'x in 
the active non-z^e fiux^. The same procedure, as applied to SNO I with SK leads to 

SNO I sinV > 0.95 at 1.23a{79% CL). (2.11) 

The upper bound on Px is therefore more restrictive (0.05) if one considers SNO I data. 
Up to 95% CL these bounds become 

sin^i; > 0.12 (SNO II) , smV > 0.36 (SNO I) (95% CL) (2.12) 

hence respectively an upper bound of 0.88 and 0.64 on the i/x fraction. These results are 
independent of the normalization to any particular standard solar model. No restriction 
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'^ Recall that the ux component is proportional to cos^ip [see eqs.(2.2), (2.3)]. 






on the sterile neutrino component, proportional to cos^a [see eqs.(2.2), (2.3)], has been 
considered so far, therefore our analysis is valid for any Vg, Vx admixture accompanying 
the LMA effect. Combining separately eqs.(2.1), (2.2) and eqs.(2.1), (2.3) one can relate 
the normalization factor Jb to the mixing angles a, '(/'■ We have respectively 

T)NC _ nCC 

Ib = i?^^ + .,,.,, j^- ^ (2.13) 

pES _ r>CC 

fB = ^^^ + . , , . 2/1- — rrv (2.14) 

We represent in fig. 3 the allowed range of sin^a (proportional to the active non- 
Ve component) in terms of fs at the 95% CL for SNO II (fig.3a) and SNO I (fig.3b) 
using inequalities (2.12) and eq.(2.13). Considering the two possible equivalent choices to 
generate fig. 3, namely eqs.(2.13) and (2.14), the former should in fact be preferred since L^ 

it leads to the narrowest error bars in Jb- The dashed and full lines in fig. 3 correspond i~n 

respectively to sin?^ = 1 (no antineutrinos P^,) and to the 95% CL upper limits for u^. ( — *v 

These are sir^il) = 0.12 for SNO II and sin?"ijj = 0.36 for SNO I [see eq.(2.12)]. Hence 
the existing shift between each adjacent dashed and full line represents the small change 
in the sterile neutrino component, proportional to cos^a, resulting from introducing in the """^ 

scheme a i>x component up to its 95% CL upper bound. This shows that the possible L\J 

sterile neutrino flux is hardly sensitive to the presence of antineutrinos, a fact whose origin (^ ^ 

becomes clear on examination of the denominator in eq.(2.13): the multiplier of sivP'a is /- — s 

very close to unity for any value of ip owing to the fact that f^ ~ 1. From fig. 3 it is also 
seen that in the absence of 9^ (x ^ e, sin'^i/j = 1) the fraction of solar neutrinos oscillating 
to active ones is greater than 0.59 (SNO II) and 0.63 (SNO I) at 2a of the non-z^g flux. 
Allowing for non-electron antineutrinos up to their 2a upper bound this fraction becomes C ") 
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J^ 



respectively 0.62 and 0.66. This result is consistent with the result of ref. [15| where the ^~^ 

authors also included KamLAND data in their analysis but were restricted to the case ^ — v 

Sin yj = 1. 

We now take an alternative view by considering separately the cases in which either 
only i^x or Vs is present along with LMA and derive in each the corresponding constraints 
on the SSM normalization factor fs- We start with the case where no steriles are present 
(only i^x)- Here siv?a = 1 and from eqs.(2.2), (2.3) one obtains 

^cc , (^^^ - R^^'W - r) - jR^'' - i^^^)(l - r-,) 

JB = a H ^7 zx TTT ^T (2.15) 

rd(r -r) - r(l - r^) 

which for SNO II and SNO I give respectively, using table I (to la) 

/b = 0.86 ±0.12 (SNO II), /ij = 0.88 ±0.13 (SNO I). (2.16) 

We note that these correspond to the allowed ranges within the lines sin?a = 1 in the two 
panels of fig. 3, the slight discrepancies with this figure being of course the result of the 
experimental uncertainties and the different procedures used for generating the two sets of 
results. For the combined SNO and SK data, eq.(2.16) becomes instead (to la) 

/b = 0.80 ±0.09 (SNO H±SK), /s = 0.84 ± 0.10 (SNO I ± SK), (2.17) 



the smaller error resulting from the smaller SK error. All these parameter ranges lie within 
the allowed la SSM error of 23% [^]. It is thus seen that the former general analysis 
which includes antineutrinos and steriles, and whose results are summarized in fig. 3, leads 
to more precise predictions for Jb, as only two experimentally measured quantities R^'-' 
and R^'^ are used in contrast to eq.(2.15). In fact, in fig. 3, where all quantities are allowed 
to vary within their 2a ranges, we have (for sin? a = 1) 

/b = 0.87 ± 0.15 (SNO II), /b = 0.91±0.19 (SNOI) (2.18) 

to be compared with eq.(2.16) where only la ranges are allowed. 

We now briefiy refer to the other special case, namely the absence of antineutrinos: v | 

only steriles are present here along with the LMA effect, hence sin^tp = 1. This case ] | ] 

corresponds to the shaded areas in fig. 3 limited by the two dashed lines and, in contrast hH 

to the previous one, no model independent equation can be obtained for /b, but only a 

degeneracy relation between /b and sirfia. This can be expressed by either of the two 

equivalent equations 

^cc ^ 

sm^a 

nES nCC 

fB = ii^^ + ^ :f (2.20) INJ 

o 



o 

nNC nCC ,'— >, 

fB = ii^^ + ^ „r,^ (2.19) Q 



r sin^a 



3. Introducing x^ Analysis 
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which correspond to eqs.(2.13) and (2.14) with sin'^'ip = 1. As previously discussed in 
the general case (LMA + i^x + t^s) the main result here is an upper bound on the sterile 

component. At 2a this is cos^a < 0.41 (from SNO II) or 0.38 (from SNO I) of the non-i^e l-t— ^ 

flux for fs = 1. ^" — ^ 

o 
o 



We refine our results by performing a x^ analysis of all four cases considered. The x^ ^ ^ 

definition is quite simple |l^ 

^,_^(ft_|^ (3.1) 

i * 

where the sum extends over the four experiments {i = ESsk,ESsnOt^C,CC), Ri, 5Ri 
denote the experimental reduced rates and their errors quoted in table I, and Rj^ are given 
by eqs. (2.1)-(2.3). The result of the x"^ minimization is shown in tables II, III for SNO II 
and SNO I respectively. 





fB 


P 

^ ee 


sin^a 


sin?ip 




LMA (2 dof) 


0.876 


0.356 


1.0 


1.0 


1.67 


LMA+P^ (1 dof) 


0.876 


0.356 


1.0 


1.0 


1.67 


LMA+z/s (1 dof) 


0.961 


0.324 


0.869 


1.0 


1.67 


LMA+z^a; + i^s (0 dof) 


0.989 


0.315 


0.833 


1.0 


1.67 





fB 


Jee 


sin^a 


sin'^tp 


Armra 


LMA (2 dof) 


0.965 


0.304 


1.0 


1.0 


2.47 


LMA+P^ (1 dof) 


0.965 


0.304 


1.0 


1.0 


2.47 


lAiK+Vs (1 dof) 


0.965 


0.304 


1.0 


1.0 


2.47 


LMA+J^^ + Us (0 dof) 


0.969 


0.302 


0.933 


1.0 


2.47 



Table III - Results of x^ analysis for SNO I. 

Finally, the 2a upper bound on the sterile component mentioned at the end of section 
2 is also shown in the contour plots of fig. 4 and corresponds to the lower edge of the light 
shaded area. 

4. Conclusions 

To conclude, the results of this paper can be summarized in figs. 2, 3 and eqs. (2.16), 
(2.17), (2.18). Allowing for antineutrinos v'x other than i>g and sterile neutrinos, both 
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Table II - Results of x^ analysis for SNO II. 

Inspection of table II (second row) shows that the best fit for case LMA+i/a; corresponds v — | 

to the very absence of Vx {sin'^ip = 1). It is also seen that allowing for Vg alone in addition | ] 

to LMA (third row) as well as LMA+i>x + i^s (fourth row) leads to a best fit solution with hj-l 
a small although non negligible i/g component (13% and 17% respectively). Furthermore 
table II also shows that Xmin ^^ independent of the values of fs, Pee, sin^a. However it 
depends on sin^ip: if in fact we let sin^ip to be unconstrained, an absolute x^in is obtained 

for an unphysical value of sin'^tp, ^. As long as sin'^ip remains constrained to its physical C_) 

region (0 < sin^ip < 1), X^m ^^ fixed regardless of the values of the other three parameters ^,_^ 

/S) Pee, sin^a. A similar situation is observed in SNO I (see table III) with the sterile |\ "v 
component totally missing {sin^a = 1) in the UAK+Vs case. This reflects the fact that the 
parameters fs, Pee, sin^a can be eliminated from eqs. (2.1) - (2.3) so as to express sin'^ip 
[vx component) in terms of experimentally measured quantities only (see eqs. (2.8), (2.9)). 

Likewise the bounds on sin'^ip are independent of the Ug component and of the other two h— i 

parameters fs. Pee (see eqs. (2.10), (2.12)). ■^— «^ 

o 
o 
o 
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^These are xLn = 0.384, sirT'ip = 2.84 (SNO II) and xLn = 0.133, sin^ip = 3.17 (SNO I). 



possibly generated in the sun through spin flavour precession accompanying LMA, we have 
derived in a model independent way, using SNO and SK data, upper bounds on the flux 
of these solar antineutrinos and steriles. We related these bounds to the parameter ranges 
allowed for the SSM normalization factor fs- To summarize: 

(i) We found an upper bound for u^ which at 2a is 0.88 (SNO II) or 0.64 (SNO I) of 
the active non-z/g flux [see fig. 2 and eq.(2.12)]. This is independent of the sterile neutrino 
component. 

(ii) In the no sterile case we obtained upper and lower bounds on fs [eqs.(2.16), (2.17), 
(2.18)]. 

(in) In the no v'x case (only steriles accompanying LMA) the fraction of solar neutrinos (~ , 

oscillating to active ones was found to be greater than 0.59 (SNO II) or 0.63 (SNO I) of t—^-i 

the non-i/e flux, a result consistent with ref. [15| which is in fact an upper bound on Vg. r~] 

(iv) Allowing, in the preceding situation, for Ux up to its 2a upper bound, these limits <- -' 

are increased by only 5%, (decrease on Us upper bound) which shows how the possible Vs vJ 

flux is hardly sensitive to the i>x component. C ^ 

(v) x^ analysis shows that the most disfavoured case (if not excluded) is Ux either with (|"~^ 

LMA or with LMA+z^j,- In SNO II it is seen that some possibility is left for LMA+z^^. 
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Figure 1: The coloured areas are the regions allowed by the la ranges of the reduced rates, as 
reported by the SNO II experiment, and sin^ip, proportional to the neutrino component of the 
active non-v^ flux. The dark shaded area is the region allowed jointly by the SNO II and SK data 
on the electron scattering reduced rate. 
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Figure 2: Same as fig.l with a magnified horizontal scale to show the lower bound on sin'^ip 
implied by the data from the two experiments. This bound is independent from the sterile neutrino 
component (see also the main text). 
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Figure 3: The allowed range at 2a of sin?' a (proportional to the total active non-i/^ flux) and Jb, 
the SSM normalization factor, using SNO II and SNO I data. Dashed lines correspond to absence 
of antineutrinos, while full lines to their upper bound at 2a. It is seen that the sterile component, 
proportional to cos^a, is hardly affected by the presence of antineutrinos. 
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Figure 4: la and 2a contour plots in the Jb, sin^a plane for SNO II data. The star represents 
the best fit point. 
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